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FIRDAT  is  a  FORTRAN  IV  program  to  compute  the  daily  components 
and  indexes  of  the  National  Fire-Danger  Rating  System.  FIRDAT  will  also 
compute  and  print  the  absolute,  relative,  and  cumulative  frequencies  of 
occurrence,  and  print  a  cumulative  frequency  distribution  for  each  of  the 
components  and  indexes. 

A  description  of  the  program  and  subroutine  along  with  examples  of 
input  and  output  data  are  presented  in  this  Note. 

Keywords:    Fire-Danger  Rating,   fire  planning,  fire  management,  data 
processing. 


Levels  for  Specific  Action  and  Manning 
Classes  for  the  new  National  Fire-Danger  Rating 
(NFDR)  System  (Deeming  et  al.  1972)  are 
determined  from  the  cumulative  frequency 
distribution  of  Burning  Index  (BI).  To  obtain 
this  distribution,  the  BI  must  be  computed  for 
all  the  available  historic  fire  weather  data.  The 
computed  BI's  then  have  to  be  sorted  into 
classes  (or  ranked)  and  the  cumulative  frequency 
computed.  This  operation  normally  would  be 
quite  time  consuming.  On  a  computer,  however, 
these  tasks  can  be  accomplished  accurately  in 
a  few  seconds  with  minimum  effort.  Using  the 
same  approach,  the  effects  of  different  land 
management  practices  or  potential  fire  situations 
are  easily  simulated,  that  is,  chipping  versus 
the  alternative  of  lopping  and  scattering  or 
directional  falling  versus  random  falling,  and 
so  forth. 

^Associate  Meteorologist,  Rocky  Mountain  Forest  and 
Range  Experiment  Station,  with  central  headquarters 
maintained  at  Fort  Collins,  in  cooperation  with  Colorado 
State  University. 

"^Systems  Analyst,  Pacific  Southwest  Forest  and  Range 
Experiment  Station,  Forest  Fire  Laboratory,  Riverside, 
California. 


This  report  documents  computer  program 
FIRDAT  constructed  to  perform  these  tasks. 
It  is  to  be  used  in  the  implementation  of  the 
NFDR  System,  and  as  an  aid  for  fire  planning 
on  the  National  Forests.  Versions  of  FIRDAT 
are  available  for  CDC  6400,  3300,  3100,  UNIVAC 
1108,  and  IBM  System/360.^ 


Program  FIRDAT 

FIRDAT  is  the  name  of  the  program  written 
in  FORTRAN  IV,  which  provides  a  fast,  accur- 
ate means  of  computing  fire  behavior  indexes 
from  historic  fire  weather  records.  It  will  pro- 
vide distributions  of  indexes  by  tabulating 
absolute,  relative,  and  cumulative  frequencies. 
A  plot  of  the  cumulative  frequency  distribution 
is  also  provided. 


^The  use  of  trade  and  company  names  is  for  the 
benefit  of  the  reader;  such  use  does  not  constitute  an 
official  endorsement  or  approval  of  any  service  or  prod- 
uct by  the  U.  S.  Department  of  Agriculture  to  the 
exclusion  of  others  that  may  be  suitable. 
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Input 

The  "data  deck"  consists  of  three  parts: 
(1)    First  Lead  Card  (see  appendix  I). — 

This  card  contains  information  about  the  ob- 
serving station,  data  input  device  and  the 
output  and  data  input  options.  The  station  in- 
formation includes  station  name,  station 
number,  slope  class,  and  elevation.  The  data 
input  device  will  designate  whether  the  ambient 
weather  data  are  available  on  punched  cards 
(unit  5)  or  magnetic  tape  (unit  1-4,  6-9).  The 
output  options  consist  of  listing  all  the  weather 
data  with  the  computed  indexes  and  compon- 
ents for  each  day;  tabulating  the  frequencies 
for  each  of  the  indexes  or  components  and 
plotting  a  graph  of  the  cumulative  frequency 
distribution;  and  punching  or  writing  on  mag- 
netic tape  the  date  and  values  of  all 
components  and  indexes.  The  components 
referred  to  are  ignition  probability,  rate  of 
spread,  and  energy  release;  the  indexes  are 
Occurrence  Index,  Burning  Index,  and  Fire 
Load  Index.  These  options  will  be  discussed 
further  in  the  Output  section. 

There  are  five  input  data  options  to  the 
program.  The  option  chosen  will  depend  upon 
the  kind  of  data  available.  (The  user,  not  the 
program,  chooses  the  option.)  The  data  com- 
binations which  determine  the  option  to  be 
used  are  tabulated  in  table  1.  Since  the  number 
of  assumptions  the  program  makes  to  satisfy 
the  data  requirements  increases  as  the  number 
of  data  parameters  decreases,  every  effort 
should  be  made  to  use  the  most  comprehensive 
data  option. 


Table  1. — Data  options  required  for  specified 
levels  of  accuracy  for  each  model — 
historical  analysis 


ua  ca 

T 

TT 
J-  -L 

options 
III  IV 

XT 
V 

Station  number 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Slope  class 

X 

X 

X 

X 

X 

uatc   \.yr  "  J  nio  .  ,   aay  j 

X 

X 

X 

X 

X 

otauc  OL  weacner 

X 

X 

X 

X 

X 

Herb,  veg .  condition 

X 

X 

X 

X 

X 

Dry— bulb  temperature 

X 

X 

X 

X 

X 

Wet-bulb  temperature 

X 

X 

X 

X 

xl 

Windspeed 

X 

X 

X 

X 

X 

Precip.  amount 

X 

X 

X 

X 

X 

24-hr.  max.  temp. 

X 

X 

X 

24-hr.  min.  temp. 

X 

X 

X 

24-hr.  max.  RH 

X 

X 

24-hr.  min.  RH 

X 

X 

1/2-inch  stick  moisture 

X 

X 

-•-If  not  available,  relative  humidity  (RH) 
or  dewpoint  may  be  substituted. 


(2)    Second  Lead  Card  (see  appendix  I).— 

This  card  contains  information  about  the  fuel 
model  and  fire  season.  Included  on  this  card 
are  the  fuel  model  designation,  loadings  in  tons 
per  acre  of  the  1-,  10-,  and  100-Hour  Timelag 
dead  fuels  and  live  fuels,  fuel  surface-area-to- 
volume  ratios  for  the  different  fuel  classes, 
the  fuel  depth,  and  the  beginning  and  ending 
dates  of  the  fire  season.  This  information  (ex- 
cept the  fire  season)  for  each  fuel  model  is 
contained  in  table  2. 


Table  2.— NFDR  fuel  models  (1972) 


Fuel  loading  Surface  area  to  volume 


Fuel 


.  model 

1-hr 

10-hr 

100-hr 

Living 

1-hr 

10-hr 

100-hr 

Living 

depth 

load 

load 

load 

load 

sigma 

sigma' 

sigma 

sigma 

Tons 

per  acre 

Feet"^ 

Feet 

A 

1.25 

0.0 

0.0 

0.0 

3000 

109 

30 

1500 

0.75 

B 

5.00 

4.0 

2.0 

2.0 

2000 

109 

30 

1500 

6.00 

C 

1.50 

1.0 

0.0 

0.0 

2700 

109 

30 

1500 

1.00 

D 

1.50 

2.5 

2.0 

0.0 

1750 

109 

30 

1500 

2.50 

E 

1.50 

1.0 

0.0 

0.0 

2500 

109 

30 

1500 

0.30 

F 

1.00 

0  .5 

0.0 

2.0 

1500 

109 

30 

1500 

2.00 

G 

3.00 

2.0 

5.0 

0.0 

1500 

109 

30 

1500 

1.25 

H 

1.00 

1.0 

1.0 

0.0 

2000 

109 

30 

1500 

0.40 

I 

4.00 

5.0 

10.0 

0.0 

1500 

109 

30 

1500 

3.50 
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(3)    Ambient  Conditions  (see  appendix  11). 

—  These  data  contain  records  of  specific  ambient 
weather  conditions  necessary  to  determine  the 
fire-danger  indexes.  Only  those  data  falling 
inside  the  fire  season  defined  on  the  Second 
Lead  Card  will  be  used  in  computing  the  fire- 
danger  indexes.  The  data  format  for  the  ambient 
conditions  (appendix  II)  was  designed  for  use 
with  the  Form  WB  612-17,  the  10-Day  Fire 
Danger  and  Weather  Record.  Depending  on  the 
input  device  designated  on  the  First  Lead  Card, 
the  program  is  capable  of  reading  the  ambient 
conditions  on  either  punched  cards  or  a  BCD 
tape  of  single  card  records.  There  is  no  limit  to 
the  number  or  records  that  can  be  read.  If 
punched  cards  are  used,  a  card  containing 
999999  in  columns  1-6  must  be  the  last  card  for 
each  station.  Computations  for  a  given  station 
will  be  completed  and  a  new  set  of  lead  cards 
will  be  read  when  the  999999  card  is  en- 
countered. If  the  data  are  on  magnetic  tape,  a 
new  set  of  lead  cards  will  be  read  when  a 
station  number  of  greater  magnitude  is  encoun- 
tered. Therefore  the  data  must  be  sorted  in 
ascending  order  by  station  number.  Fuel  mois- 
ture values  carry  over  from  day  to  day,  so  the 
data  should  be  arranged  in  ascending  order  by 
date  for  each  station  on  both  cards  and  tape. 
See  figure  1  for  an  example  of  an  input  data 
deck. 


Output 

There  are  four  parts  to  the  output  of  pro- 
gram FIRDAT: 

(1)  Daily  Information  List. —  This  part  of 
the  output  provides  a  line  per  card  listing  of 
weather  parameters  and  fire-danger  components 
and  indexes.  It  is  optional  and  may  be 
suppressed  by  making  the  proper  entry  on  the 
First  Lead  Card. 

(2)  Frequency  Tabulations.  —  For  any  or 
all  of  the  six  fire-danger  components  and 
indexes,  the  absolute,  relative,  and  cumulative 
frequencies  are  tabulated  along  with  the  class 
upper  boundary.  The  choice  of  indexes  to  be 
tabulated  is  made  on  the  First  Lead  Card  (see 
appendix  I). 

(3)  Cumulative  Frequency  Distribution. — 
A  cumulative  frequency  curve  will  be  printed 
on  the  computer  page  for  the  indexes  designated 
in  the  output  options. 

(4)  Auxiliary  Output  of  Components  and 
Indexes.  — The  program  has  an  output  feature 
whereby  the  station  number,  date,  and  the  six 
fire  components  and  indexes  can  be  punched  on 
cards  or  written  in  BCD  on  magnetic  tape.  This 
option  is  on  the  First  Lead  Card. 


The  auxiliary  output  will  contain  two  header 
cards  per  station  run  plus  one  card  per  observa- 
tion containing  indices,  arranged  and  formated 
as  follows: 


Header  Card  1  Col.  1-6  Station  Number 
Header  Card  2    Col.     1-6     Station  Number 

7-10  Fuel  Model  Type 
1  1-30  Station  Name 

3  1  Slope  Class 
32-36  Station  Elevation  (feet) 
Card/observation    1-6     Station  Number 
7-8  Year 
9-10  Month 
11-12  Day 

13-15  Probability  of  Ignition 
16-18  Spread  Factor 
1  9-2  1   Energy  Release 
22-24  Occurrence  Index 
25-27  Burning  Index 
28-30  Fire  Load  Index 


This  output  is  arranged  to  be  compatible 
with  a  foUowup  program  which  will  display 
the  seasonal  variation  of  the  indexes  graphically. 

Examples  of  the  output  are  found  in 
appendix  III  ;  a  flow  chart  of  FIRDAT  is  in 
appendix  IV. 

Assumptions 

For  many  initial  applications  of  the  NFDR 
System,  the  historic  data  available  will  not  in- 
clude some  of  the  parameters  needed  to  compute 
the  fire-danger  indexes  and  components.  Some 
assumptions  were  made  to  fulfill  data  require- 
ments. These  assumptions  will  depend  upon 
the  data  available  which,  in  turn,  will  determine 
which  data  option  will  be  chosen  by  the  user. 
The  following  assumptions  were  made  for  each 
option. 

Option  I.  —  No  assumptions  are  made  to  fill 
missing  data.  All  necessary  data  are  provided. 

Option  II.  — The  10-Hour  Timelag  Fuel  Mois- 
ture is  computed  from  observed  temperature 
and  humidity  instead  of  being  observed  from 
1/2-inch  fuel  moisture  sticks  (analog). 

Option  III.  —  In  addition  to  the  assumption 
in  Option  II,  the  maximum  and  minimum  rela- 
tive humidities  are  computed  from  the  maximum 
and  minimum  temperatures,  assuming  constant 
specific  humidity  for  the  day. 
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Option  IV.  — The  100-Hour  Timelag  Fuel 
Moisture  is  computed  from  "yesterday's"  10- 
and  100-Hour  Timelag  Fuel  Moisture  and 
''today's"  10-Hour  Timelag  Fuel  Moisture  rather 
than  from  average  temperature  and  humidity. 
The  10-Hour  Timelag  Fuel  Moisture  is  the 
measured  fuel  moisture  from  the  1/2-inch  fuel 
moisture  sticks. 

Option  v.  — This  option  is  a  combination  of 
Options   n    and  IV. 

In  addition  to  the  assumptions  made  as 
required  by  the  options,  other  assumptions  are 
made  for  each  observation  and  are  independent 
of  which  option  is  selected.  Among  these  are: 


a.  All  risk  values  are  assigned  inside  the  pro- 
gram. Lightning  and  man-caused  risk  are 
set  to  5  each  for  a  total  risk  of  10. 

b.  If  precipitation  has  occurred  in  any  amount 
during  the  past  24  hours,  the  precipitation 
duration  is  set  to  1  hour. 

c.  If  the  state  of  the  weather  indicates  precipita- 
tion was  occurring  at  the  time  of  the  obser- 
vation, the  fine  fuel  moisture  and  the  10-Hour 
Timelag  Fuel  Moisture  (except  for  1/2-inch 
fuel  sticks)  is  set  to  31  percent  and  25  per- 
cent, respectively. 

d.  Living  fuel  moisture  (woody)  is  set  to  100 
percent. 

In  selecting  these  options  we  felt  that  the 
five  most  desirable  and  likely  combinations  of 
data  parameters  were  represented. 


Missing  Data 

At  several  points  in  the  program,  certain 
data  parameters  are  checked  to  determine  if 
they  have  a  zero  value  or  are  missing  (blank). 
If  zero  or  missing  values  are  detected,  that  data 
card  is  usually  ignored  and  the  program  pro- 
ceeds to  the  next  card.  The  parameters  checked 
include: 


Value  to 

Parameter 

reject 

Option 

Temperature 

m  isg/u'' 

l-V 

Relative  humidity 

0/misg/u 

l-V 

Windspeed 

misg/u 

l-V 

Tmax  -  Tmin 

misg/u 

l-lll 

RHmax  -  RHmin 

m  isg/u 

1,11 

1  /2-inch  stick  moisture 

0/misg 

l,IV 

"u  -  unrecognizable  character. 


Subroutines 

Unless  otherwise  stated,  all  subroutines  and 
functions  are  called  only  from  the  main  program. 
Subroutines  are  arranged  in  alphabetical  order, 
not  in  the  order  in  which  they  are  called. 


Function  ADJFFM 

Function  ADJFFM  will  return  a  value  for 
fine  fuel  moisture  in  percent  for  values  of 
1-Hour  Timelag  Fuel  Moisture  and  percent  green 
herbaceous  vegetation.  This  function  is 
essentially  a  table  look-up  process.  No  adjust- 
ment will  be  made  to  the  1-Hour  Timelag  Fuel 
Moisture  if  the  herbaceous  vegetative  condition 
is  cured  (less  than  10  percent  green  material, 
appendix  II). 

An  unrecognizable  value  for  herbaceous 
condition  or  stage  will  result  in  a  default  value 
of  30  percent  for  percentage  of  green  herbaceous 
material. 


Subroutine  C0DPRC 

This  routine  accepts  encoded  over-punched 
numeric  data,  originally  read  in  Al  format,  and 
returns  the  numeric  value  signed  plus.  A 
missing-flag  is  returned  ".TRUE."  if  the  field 
contains  any  invalid  code  combination  (legiti- 
mate codes  are  0,  1,  2,  3,  4,  5,  6,  7,  8,  9,  11-0, 

11-  1, 11-2,  etc.,  11-9,  12-0,  12-1,  .  .  .  12-0,  12-1,  .  .  . 

12-  9,  12-ONLY,  11-ONLY,  BLANK).  A  blank-flag 
is  set  ".TRUE."  if  the  field  is  blank,  and  a 
negative-flag  is  set  ".TRUE."  if  anll-overpunch 
is  present  (12-overpunch  treated  as  if  not 
present,  i.e.,  12-3  =  3).  The  returned  floating 
point  value  is  always  signed  plus.  C0DPRC  is 
called  from  subroutine  SCAN. 


Subroutine  CUM 

When  all  the  data  for  a  given  station  have 
been  read,  the  subroutine  CUM  is  called.  Here 
the  absolute,  relative,  and  cumulative  frequen- 
cies are  computed  for  the  desired  indexes,  and 
the  number  of  classes  is  decided  upon  for 
plotting  purposes.  The  number  of  classes  chosen 
is  that  necessary  to  account  for  99.5  percent  of 
the  data. 


Function  EQMC 

This  function  returns  a  value  of  equilibrium 
moisture  content  for  wood,  given  temperature 
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(°F)  and  relative  humidity  (%).  The  value  is 
derived  from  regressions  of  the  Forest  Products 
Laboratory  Tables  (Simard  1968),  This  function 
is  called  from  Subroutine  MlOO  as  well  as  from 
the  main  program. 


Function  FFMCHS 

This  function  computes  the  1-Hour  Timelag 
Fuel  Moisture  (percent).  Function  arguments 
are  dry  bulb  temperature  (°F),  relative  humidity 
(%),  and  cloud  cover  (decimal). 


Function  GFNC 

Given  the  packing  ratio  (B),  volume  to  area 
ratio  (S),  and  optimum  packing  ratio  (B0PT), 
this  function  will  deliver  the  maximum  reaction 
intensity. 


Function  H2PCNV 

This  function  returns  approximate  pressures 
in  millibars  given  elevation  in  feet  m.s.l.  The 
function  uses  the  same  elevation  ranges  used 
in  Weather  Bureau  psychrometric  tables. 


Real  Function  MXFNC 

Given  percent  green  and  1-Hour  Timelag 
Fuel  Moisture,  this  function  computes  and 
returns  the  extinction  moisture  content  (%) 
for  living  herbaceous  fuels. 


Subroutine  MlOO 

This  subroutine  computes  the  100-Hour 
Timelag  Fuel  Moisture  (%)  from  Fosberg's 
(1972)  theoretical  solution.  The  arguments  for 
this  subroutine  are  24-hour  average  temperature 
(°F),  average  relative  humidity  (%),  and  dura- 
tion of  precipitation  (hrs)  for  the  preceding  24 
hours.  This  subroutine  is  used  in  options  I-HI. 
Function  EQMC  is  called  to  deliver  equilibrium 
moisture  content. 


Subroutine  MIOOA 

This  subroutine  estimates  the  100-Hour 
Timelag  Fuel  Moisture  from  a  stochastic  model 
given  inputs  of  yesterday's  10-  and  100-Hour 
Timelag  Fuel  Moisture  and  today's  10-Hour 
Timelag  Fuel  Moisture.  This  subroutine  is  used 


in  options  IV  and  V  where  average  temperature 
and  humidity  are  not  available. 


Subroutine  PL0TR 

This  subroutine  contains  the  instructions  to 
plot  the  cumulative  frequency  distribution  of 
the  index  specified  by  the  index  number  in  the 
argument  list.  If  the  number  of  classes  is  less 
than  or  equal  to  five,  or  if  there  are  fewer  than 
100  data  points,  no  graph  will  be  plotted  and 
the  appropriate  error  message  will  be  printed. 


Subroutine  PRINTR 

Class  frequencies  are  printed  in  tabular  form 
by  this  subroutine.  The  variables  to  be  printed 
are  in  blank  "C0MM0N".  The  arguments  nec- 
essary are  the  number  which  designates  which 
index  is  being  handled  and  the  number  of 
classes  to  be  considered. 


Subroutine  PSYCHR 

This  subroutine  computes  relative  humidity 
(%)  from  dry  bulb  temperature  (°A),  wet  bulb 
temperature  (°A),  and  pressure  (mb).  Sub- 
routine VAP0R  is  called  from  PSYCHR  to 
deliver  saturation  vapor  pressures.  Wet  bulb 
temperatures  below  0°  C.  may  be  used  but  no 
correction  has  been  made  for  wet  bulb  icing. 
The  humidities  returned  will  be  between  1  and 
100  percent  inclusive. 


Subroutine  SCAN 

This  routine  performs  alpha-string  to  float- 
ing-point conversion.  "-ALPHA-"  is  an  alpha- 
numeric array  read  as  N-Al,  "-VALUE-"  is  the 
returned  real  value,  "-NWID-"  is  the  number 
of  characters  of  the  field  width,  "-NDEC-"  is 
the  number  of  digits  to  be  considered  to  be  to 
the  right  of  the  decimal  point  (same  con- 
ventions as  FORTRAN  format  specifications), 
"-MISFLG-"  is  a  logical  set  ".TRUE."  for 
invalid  code  combinations  or  illegal  sequences, 
"-BLKFLG-"  is  a  logical  set  ".TRUE."  if  the 
field  in  question  is  entirely  blank.  "-0VP0S-" 
is  an  integer  specifying  which  character  of  the 
field  may  be  legally  overpunched.  For  instance, 
for  a  4  character  field,  if  "-OVPOS-"  were  3, 
then  if  an  overpunch  (11-PUNCH)  were  found 
in  any  position  except  the  3rd  character,  the 
error-flag  would  be  set  ".TRUE."  and  the  field 
considered  invalid.  A  "0"  value  for  "-0VP0S-" 
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means  that  any  character  may  be  overpunched, 
while  any  negative  value  of  "-OVPOS-"  means 
that  no  character  may  be  overpunched. 
"-MINFLG-"  is  set  ".TRUE."  if  the  field  is 
negative  (overpunched).  "-MINFLG-"  is 
returned  to  permit  the  situation  of  an  over- 
punched  ''0"  field.  CODPRC  is  called  from  SCAN. 


Subroutine  SET 

In  "subroutine  SET,  the  arrays  in  which  the 
frequencies  are  accumulated  are  initialized  and 
the  class  interval  is  determined.  This  subroutine 
is  used  only  once  at  the  beginning  of  each 
station. 


Subroutine  SORT 

In  subroutine  SORT,  the  index  value  is 
sorted  into  its  proper  class,  and  the  number-in- 
class  counters  are  incremented.  This  subroutine 
is  accessed  for  every  data  card. 


Function  U20FNC 

This  function  reduces  wind  measured  at 
20  feet  aboveground  to  midflame  height.  It  also 
restricts  the  maximum  reduced  wind  to  limits 
required  by  theory  for  the  fuel  bed. 


Subroutine  VAP0R 

For  a  temperature  ("A),  this  subroutine 
will  return  the  saturation  vapor  pressure  over 
water  (mb).  The  operation  is  based  on  the 
equation  from  List's  (1966)  Smithsonian 
Meteorological  Tables  (p.  350).  This  subroutine 
is  called  from  program  FIRDAT  and  subroutine 
PSYCHR. 
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Appendix  I:   Lead  card  formats 


Card  field 


Columns 
From  To 


Number 
of 

columns 


Remarks 


First  Lead  Card  -  Station  parameters 


Station  name 
Station  number 
Slope  class 
Elevation 
Input  device 
designation 


Indexes  output  device 
designation 

Ignition  probability 
Rate  of  spread 
Energy  release 
Occurrence  index 
Burning  index 
Fire  load  index 
Daily  info  list 
Dota  option 


21 


34 


20 
26 
31 
38 

41 


42 

44 
46 
48 
50 
52 
54 
56 
60 


20 
6 


Feet  (right  justify) 

1  -4,  6-9  =  Magnetic  tope 
input  (assign  magnetic 
tape  logical  unit) 

5  =  Cord  reader  input 

0  -  No  output 

1-9  Fortran  logical  file  no. 

Enter  T  to  suppress^ 


-5  Inclusive  (see  table 


Second  Lead  Cord  -  Model  name  and  fuel  bed  parameters 


Model  type 

Fuel  loadings 

1  -hour  timelag 
1  0-hour  timelag 
1  00-hour  timelag 
Live 

Area  volume  ratio 
1  -hour  timelag 
1  0-hour  timelag 
1  00-hour  timelag 
Live 


5 
10 
15 
20 

25 
30 
35 
40 


9 
14 
19 
24 

29 
34 
39 
44 


(See  table  2) 

Tons  per  acre 

Right  justify  in  field' 


Bed  depth 

Fire  season 

Beginning  month 
Beginning  day 
Ending  month 
Ending  day 


45 


50 
52 
55 

57 


49 


51 
58 
56 
58 


Integer. 

For  CDC  3100  version,  use  1  instead  of  T. 
'Real  number  —  right  justify  or  punch  decimal. 
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Appendix  II:   Card  punching  and  verifying  instructions  for  the  weather  records 


Number 


Field 

Card  field 

Columns 

of 

Remarks 

no. 

From 

To 

columns 

Stotion  number 

1 

6 

6 

upper  iiyHi  corner 

Hour  (observation  time) 

7 

8 

2 

First  two  digits 

Year  (period  of  time) 

9 

10 

2 

Last  two  digits 

Month  (period  of  time) 

1  1 

12 

2 

(Convert  A  months  to  N) 
(Jan  01;  Feb  02,  etc.) 

1 

Day 

1  3 

14 

2 

2 

State  of  weather 

15 

15 

1 

3 

Herbaceous  stage 

16 

16 

1 

(M 

4 

Dry  bulb  temp. 

17 

19 

3 

Minus,  "-"  overpunch  col. 

19 

5 

Wet  bulb  temp. 

20 

21 

2 

Minus,       overpunch  col. 

21 

7 

Dew  point  temp. 

24 

25 

2 

Minus,  "-"  overpunch  col. 

25 

8 

Relative  humidity 

26 

27 

2 

100%,  punch  zeros,       over  26 

14 

Windspeed 

35 

36 

2 

100  &  over,  punch  zeros, 
over  35 

'/a"  fuel  moisture  sticks 

44 

46 

3 

1  decimal  place 

34 

Maximum  temp. 

59 

61 

3 

Minus,  "-"  overpunch  col. 

61 

35 

Minimum  temp. 

62 

63 

2 

Minus,       overpunch  col. 

63 

36 

Maximum  humidity 

64 

65 

2 

100%,  punch  zeros  "-"  over  64 

37 

Minimum  humidity 

66 

67 

2 

41 

24-hour  amount 

73 

75 

3 

2  decimal  places.  Decima 

blank,  punch  zero.  Trace,  punch 


TOO 


No  data  in  field,  leave  blank. 
Left  fill  all  fields  with  zero. 

Punch  Operator  Functions: 

The  field  numbers  on  the  card  punch  instructions  correspond  with  the  numbers  printed  on  the  WB 
Form612-17. 

One  detail  card  per  day  of  the  month  entry  will  be  punched  and  verified. 

Columns  1  through  12  can  be  duplicated  for  each  form  after  the  first  day  is  punched. 

Credits  or  minus  figures  will  be  punched  with  a  "-"  overpunch  in  the  units  (low  order)  position. 

Figures  which  exceed  two  digits,  100  or  more,  will  be  punched  with  a  "-"  overpunch  in  the  high 
order  position. 

Fields  without  data,  leave  blank. 
Left  fill  all  data  fields  with  zeros. 


1-cured  (10  percent  green  material);  2-transition  (30  percent  green  material);  3-green  (50  percent 
green  material). 
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Appendix  III:  Output 

Abbreviations  used  on 
Daily  Information  List 


SW  -  state  of  the  weather  L  F 

D  B  TMP  -  dry  bulb  temperature  PC  GR 

REL  HUM  -  relative  humidity  1-HR  FM 

WD  SP  -  windspeed  ADJ  1-HR 

TMP  MAX  -  24-hour  maximum  temperature  10-HR  FM 

TMP  MIN  -  24-hour  minimum  temperature  100-HR  FM 

R  H  MAX  -  24-hour  maximum  relative  humidity  IGN  COM 

R  H  MIN  -  24-hour  minimum  relative  humidity  SPD  COM 

PPT  AMT  -  24-hour  precipitation  amount  ENRL  COM 

PPT  DUR  -  precipitation  duration  (either  0  or  1 )  OCC  NDX 

L  R  -  Lightning  Risk  BUR  NDX 

MR  -  Man-Caused  Risk 


living  fuel  moisture  (woody  vegetation  condition) 

percentage  green  herbaceous  material 

1-hour  timelag  dead  fuel  moisture 

fine  fuel  moisture 

10-hour  timelag  dead  fuel  moisture 

100-hour  timelag  dead  fuel  moisture 

Ignition  Component 

Spread  Component 

Energy  Release  Component 

Occurrence  Index 

Burning  Index 


NATIONS,  rise  CANCER  RATINC  STSTEH  -  OPTION  1  05/21/72 
STATIJM  -  ciL*  CENTER  (2920061  TuCL  MOOEU  C  SLOPE    3        EtEV.  S700 


U  A  ■  L 

c 

9 

V  D 

DCl 

•A/ 

R  H 

R  H 

PPT 

PPT 

L     M  L 

PC 

1  -VR 

AOJ 

1 

1 

3^ 

OCC 

w 

mjM 

MAX 

MIN 

MAX 

MIN 

AMT 

R  R  r 

CR 

TM 

1 

FM 

COn 

COM 

»*x 

NDX 

71 0418 

5 

59 

92 

€ 

69 

57 

100 

31 

010 

0. 

5    5  7 

10 

31 .0 

50.0 

29.0 

0. 

0 

0 

0 

0 

0 

7  I  v«  1  y 

«^ 

3 ' 

m* 

•  • 

44 

27 

100 

62 

012 

0. 

5    9  7 

10 

19.0 

19.0 

15.9 

0 . 

0 

9 

6 

0 

6 

71 0^0 

0 

54 

92 

1 4 

54 

19 

100 

35 

001 

0. 

9    5  7 

10 

6.0 

6.0 

7.9 

0. 

46 

1 1 

20 

47 

14 

710^  1 

5 

50 

92 

1  5 

99 

21 

100 

23 

000 

0. 

9    9  7 

10 

7.0 

7.0 

9. 0 

0. 

55 

9 

21 

56 

15 

A 

0 

9V 

<4 

1  4 

60 

29 

77 

23 

000 

0. 

1 0 

4.0 

4,0 

9.5 

0 . 

64 

1 4 

26 

66 

1 7 

710*23 

0 

67 

24 

1 1 

66 

23 

100 

23 

000 

0. 

5    5  7 

to 

4.0 

4.0 

2.0 

0. 

66 

9 

29 

99 

19 

710424 

c 

67 

20 

1 1 

70 

26 

too 

19 

000 

0. 

to 

4.0 

4.0 

2.0 

0. 

66 

9 

29 

69 

15 

710425 

c 

67 

17 

i« 

71 

40 

45 

19 

000 

o; 

10 

3.0 

5.0 

1.9 

0. 

76 

24 

50 

79 

25 

710426 

0 

59 

t« 

19 

70 

43 

54 

12 

000 

0. 

10 

4.0 

4.0 

1.9 

0. 

64 

21 

29 

66 

21 

710427 

0 

65 

15 

1 1 

66 

22 

74 

13 

000 

0. 

10 

3.0 

5.0 

1.0 

0. 

76 

10 

50 

79 

16 

71042« 

0 

63 

19 

13 

69 

24 

99 

16 

000 

0. 

10 

4.0 

4.0 

1.0 

0. 

65 

12 

29 

67 

17 

71042* 

AN 

ERROR 

HAS 

BEEN 

rouND 

IN 

T>C  DATA. 

RECORD 

HAS 

REJECTED. 

710430 

AN 

ERROR 

HAS 

BEEN 

rouND 

IN 

THE  DATA. 

RECORO 

WAS 

REJECTED. 

710501 

0 

7« 

10 

9 

79 

25 

95 

10 

000 

0. 

10 

2.0 

2.0 

1.0 

0. 

01 

9 

52 

95 

15 

710502 

AN 

ERROR 

HAS 

BEEN 

rouND 

IN 

J>€.  DATA. 

RECORO 

MAS 

REJECTED. 

710503 

0 

79 

10 

15 

92 

25 

95 

9 

000 

0. 

to 

2.0 

2.0 

5.0 

0. 

91 

16 

30 

95 

19 

710504 

0 

69 

• 

19 

•6 

55 

29 

6 

000 

0. 

to 

2.0 

2.0 

5.0 

0. 

99 

27 

50 

91 

24 

710505 

0 

66 

15 

15 

71 

26 

60 

12 

000 

0. 

10 

5.0 

5.0 

5.0 

0. 

76 

14 

29 

79 

19 

710506 

2 

62 

21 

9 

69 

52 

66 

19 

000 

0. 

10 

9.0 

9.0 

5.9 

0. 

51 

6 

29 

52 

12 

710507 

2 

67 

20 

15 

70 

54 

71 

10 

000 

0. 

to 

9.0 

9.0 

4.0 

0. 

52 

1 1 

25 

55 

15 

71050* 

1 

69 

16 

16 

73 

50 

75 

15 

000 

0. 

to 

5.0 

5.0 

4.0 

0. 

77 

19 

27 

79 

20 

71050« 

0 

70 

12 

15 

71 

29 

57 

12 

000 

0. 

10 

2.0 

2.0 

5.9 

0. 

99 

20 

29 

9t 

21 

710510 

1 

70 

15 

6 

72 

27 

75 

14 

000 

0. 

to 

5.0 

5.0 

5.9 

0. 

77 

4 

29 

79 

10 

71051 1 

9 

65 

25 

15 

75 

50 

92 

15 

000 

0. 

to 

9.0 

9.0 

4.9 

0. 

92 

14 

24 

55 

17 

710512 

5 

64 

56 

7 

65 

34 

100 

52 

016 

0. 

to 

9.0 

9.0 

12.0 

0. 

52 

5 

15 

55 

7 

710513 

1 

71 

16 

5 

71 

51 

71 

16 

019 

0. 

10 

5.0 

5.0 

12.9 

0. 

77 

3 

20 

79 

9 

710514 

1 

76 

26 

7 

76 

56 

100 

16 

000 

0. 

10 

4.0 

4.0 

9.0 

0. 

69 

4 

22 

70 

10 

710515 

0 

79 

16 

1 1 

•0 

54 

96 

16 

000 

0. 

to 

5.0 

5.0 

6.9 

0. 

79 

to 

25 

91 

15 

710516 

0 

79 

15 

1 1 

•2 

55 

76 

15 

000 

0. 

10 

5.0 

5.0 

9.9 

0. 

79 

10 

26 

91 

15 

710517 

1 

71 

16 

16 

•  1 

55 

67 

12 

000 

0. 

to 

5.0 

5.0 

9.0 

0. 

77 

19 

26 

79 

20 

7105H 

0 

65 

12 

7 

74 

59 

55 

9 

000 

0. 

to 

2.0 

2.0 

9.0 

0. 

97 

6 

29 

99 

12 

710519 

0 

71 

10 

• 

71 

24 

62 

1 1 

000 

0. 

10 

2.0 

2.0 

4.0 

0. 

99 

7 

29 

91 

15 

710520 

0 

76 

<0 

1  1 

79 

24 

77 

1 1 

000 

0. 

10 

2.0 

2.0 

5.9 

0. 

91 

12 

29 

95 

17 

710921 

5 

75 

12 

25 

•0 

29 

66 

10 

000 

0. 

10 

5.0 

5.0 

5.9 

0. 

75 

37 

29 

75 

27 

710922 

0 

70 

15 

16 

79 

29 

69 

1 1 

000 

0. 

to 

5.0 

3.0 

5.9 

0. 

77 

20 

29 

79 

20 

710925 

0 

69 

16 

16 

75 

29 

75 

IS 

000 

0. 

to 

3.0 

3.0 

5.9 

0. 

77 

20 

29 

79 

20 

710924 

0 

75 

7 

4 

79 

26 

46 

7 

000 

0. 

to 

t.O 

t.O 

5.9 

0. 

too 

3 

31 

too 

10 

710925 

0 

•  1 

1  1 

9 

•  1 

29 

77 

7 

000 

0. 

10 

2.0 

2.0 

5.0 

0. 

92 

9 

50 

94 

15 

710926 

0 

•0 

10 

1 1 

•2 

52 

61 

10 

000 

0. 

to 

2.0 

2.0 

5.0 

0. 

92 

12 

50 

94 

17 

710927 

0 

•0 

12 

12 

95 

52 

72 

10 

000 

0. 

to 

2.0 

2.0 

5.0 

0. 

92 

14 

50 

94 

19 

710929 

1 

77 

14 

1 1 

92 

45 

46 

12 

000 

0. 

to 

3.0 

3.0 

5.9 

0. 

79 

10 

29 

91 

15 

710929 

2 

65 

21 

17 

79 

54 

69 

14 

000 

0. 

10 

9.0 

9.0 

5.9 

0. 

92 

17 

25 

55 

19 

710950 

0 

69 

19 

12 

70 

51 

90 

19 

000 

0. 

to 

4.0 

4.0 

5.9 

0. 

66 

1  1 

26 

69 

15 

710951 

0 

75 

15 

15 

74 

52 

66 

15 

000 

0. 

10 

3.0 

3.0 

4.0 

0. 

79 

17 

27 

90 

19 

710601 

0 

77 

9 

12 

77 

29 

57 

10 

000 

0. 

10 

2.0 

2.0 

5.5 

0. 

91 

14 

29 

95 

19 

710602 

5 

7i 

• 

14 

90 

54 

42 

9 

000 

0. 

to 

2.0 

2.0 

5.5 

0. 

96 

19 

29 

99 

20 

710605 

I 

76 

10 

It 

7* 

51 

9* 

9 

000 

0. 

to 

2.0 

2.0 

5.5 

0. 

91 

27 

29 

95 

24 

710604 

0 

77 

7 

1 1 

79 

22 

57 

7 

000 

0. 

to 

1.0 

1.0 

5.0 

0. 

100 

14 

51 

100 

1* 

710605 

0 

•0 

10 

19 

92 

24 

94 

7 

000 

0. 

10 

2.0 

2.0 

5.0 

0. 

92 

30 

50 

94 

26 

710606 

0 

•2 

• 

15 

94 

29 

96 

9 

000 

0. 

to 

1.0 

t.O 

5.0 

0. 

too 

19 

51 

100 

21 
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NATIOMAL  riRt:  DAMCER  R4TINC  SYSTtM  -  OPTION  t 


rREQUEfJCT  DISTRIBUTIONS  TOR  BURMIMC  IKIEx 


05/2' 


STATION!  -  GIL*  CENTER  -292006     *     SLOPE  CL*SS  3     '     HOOEL  C  *  SEASON    401-1051     '     NO.  PNTS.  220 

CLASS  UPPER  NO.    IN  RELATIVE  CUMULATIVE 


3. 

SOT 

CLASS 

rREQ  (PCT) 

rREQ 

1 

2 

8.6 

8.6 

2 

4 

3e 

17.5 

25.9 

3 

£ 

26 

12.7 

38.6 

i 

e 

20 

e.  1 

47.7 

5 

10 

29 

15.2 

60.9 

6 

»2 

25 

10.5 

71  .4 

7 

t4 

1 1 

5.0 

76.4 

e 

16 

21 

0.5 

65.9 

e 

ie 

15 

6,6 

02.7 

10 

20 

e 

5.6 

06.4 

1 1 

22 

5 

t.4 

97.7 

12 

24 

5 

1.4 

09.  1 

13 

26 

1 

.5 

99.5 

14 

26 

1 

0.0 

100.0 

15 

50 

0 

0.0 

100.0 

16 

52 

0 

0.0 

100.0 

17 

54 

0 

0.0 

100.0 

18 

56 

9 

0.0 

100.0 

10 

5« 

0 

0.0 

100.0 

20 

40 

0 

0.0 

100.0 

21 

42 

0 

0.0 

100.0 

22 

44 

0 

0.0 

100.0 

23 

46 

0 

0.0 

100.0 

24 

48 

0 

0.0 

100.0 

25 

50 

0 

0.0 

100.0 

NATIONAL  f  IRE  DANCER  RATINC  SYSTEM  -  OPTION  1  CUMULATIVE  TREQUEfJCT    (PCT)   -  BURNINC  INDEX 

STATION  -  GiLA  CENTER  -292006    '     SLOPE  CLASS  3    *    MODEL  C         *    f\t>Z  SEASON    401-1031     *    NO.  PtiTS.  220 

50 


40 


30 


10 


0  !-  2:  35  40  50  60  70  80  91 

PERCEN' 


1  1 


Appendix  IV:   Flow  chart  for  program  FIRDAT 


READ  LEAD 
CARDS 


COMPUTE  WEIGHTING 
FACTORS  FROM 
LOADINGS  &  VOL.- 
AREA  RATIOS 


READ  FIRE 
WEATHER  DATA 


 xz  

CONVERT 
CHARACTERS  TO 
NUMBERS 


COMPUTE 
NEEDED 

PHYSICAL  INFO. 


SL. 


COMPUTE 
FUEL 

MOISTURES 


V  

COMPUTE 
COMPONENTS  & 
INDEXES 


OUTPUT: 

DAILY  INFO.  LIST, 
PUNCH  CARDS 


Agriculture— CSU,  Ft.  Collins 


SL 


OUTPUT: 

TABULATION, 

GRAPH 


1  2 


